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Abstract 

The Weyl correspondence and the related Wigner formalism lie at 
the core of traditional quantum mechanics. We discuss here an alter- 
native quantization scheme, whose idea goes back to Born and Jor- 
dan, and which has recently been revived in another context, namely 
time-frequency analysis. We show that in particular the uncertainty 
principle does not enjoy full symplectic covariance properties in the 
Born and Jordan scheme, as opposed to what happens in the Weyl 
quantization. 

1 Introduction 

The problem of "quantization" of an "observable" harks back to the early 
days of quantum theory; mathematically speaking, and to use a modern 
language, it is the problem of assigning to a symbol a pseudo-differential 
operator in a way which is consistent with certain requirements (symme- 
tries under a group of transformations, positivity, etc.). Two of the most 
popular quantization schemes are the Kohn-Nirenberg and Weyl correspon- 
dences. The first is widely used in the theory of partial differential equations 
and in time- frequency analysis (mainly for numerical reasons), the second 
is the traditional quantization used in quantum mechanics. Both are actu- 
ally particular cases of Shubin's pseudo-differential calculus, where one can 
associate to a given symbol a an infinite family (A T ) T of pseudo-differential 
operators parametrized by a real number r, the cases r = 1 and r = | cor- 
responding to, respectively, Kohn-Nirenberg and Weyl operators. It turns 
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out that each of Shubin's r-operators can be alternatively defined in terms 
of a generalization Wig r of the usual Wigner distribution by the formula 

<A r V^) = (a,Wig r (V^)) 

and this observation has recently been used by researchers in time-frequency 
analysis to obtain more realistic phase-space distributions (more about this 
in the discussion at the end of the paper). They actually went one step 
further by introducing a new distribution by averaging Wig r for the values 
of r in the interval [0, 1]. This leads, via the analogue of the formula above 
to a third class of pseudo-differential operators, corresponding to the aver- 
aging of Shubin's operators A T . It was noted by the present author that 
this averaged pseudo-differential calculus is actually an extension of one of 
the first quantization schemes discovered by Born, Jordan, and Heisenberg 
around 1927, prior to that of Weyl's. 

The aim of the present paper is to give a detailed comparative study 
of the Weyl and Born- Jordan correspondences (or "quantization schemes" ) 
with an emphasis on the symplectic covariance properties of the associated 
uncertainty principles. 

Notation 1 We write x = (x\,...,x n ) and p = (pi,..-,p n ) a- n d z = (%,p)- 
In matrix calculations x,p,z are viewed as column vectors. The phase space 
M 2n = R n x W 1 is equipped with the standard symplectic form a(z,z') = 

px' — p'x; equivalently a(z,z') = Jz ■ z' where J = ( ^nx™ ^ xn J i s the 

standard symplectic matrix. We denote by S(M? n ) the Schwartz space of 
rapidly decreasing smooth functions and by S'(M. 2n ) its dual (the tempered 
distributions). 

2 Discussion of Quantization 

After Werner Heisenberg's seminal 1925 paper [23J which gave rigorous bases 
to the newly born "quantum mechanics" , Born and Jordan [4 J wrote the first 
comprehensive exposition on matrix mechanics, followed by an article with 
Heisenberg himself |5j. These articles were an attempt to solve an ordering 
problem: assume that some quantization process associated to the canonical 
variables x (position) and p (momentum) two operators X and P satisfying 
the canonical commutation rule XP — PX = ih. What should then the 
operator associated to the monomial x m p n be? Born and Jordan's answer 
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was 
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which immediately leads to the "symmetrized" operator ^(XP + PX) when 
the product is xp. In fact Weyl and Born- Jordan quantization lead to the 
same operators for all powers x' m or p n , or for the product xp (for a detailed 
analysis of Born and Jordan's derivation see Fedak and Prentis [9], also 
Castellani [6j and Crehan [7]). Approximately at the same time Hermann 
Weyl had started to develop his ideas of how to quantize the observables 
of a physical system, and communicated them to Max Born and Pascual 
Jordan (see Scholz [28]). His basic ideas of a group theoretical approach 
were published two years later |33|, I34j . One very interesting novelty in 
Weyl's approach was that he proposed to associate to an observable of a 
physical system what we would call today a Fourier integral operator. In 
fact, writing the observable as an inverse Fourier transform 

a(x,p)= f e i{ps+xt) Ta(s,t)dsdt (2) 

he defined its operator analogue by 

A= f e i{ps+jlt) Ta{s,t)dsdt (3) 

which is essentially the modern definition that will be given below (formula 

([9]). We will denote the Weyl correspondence by a <— ^ A^j or Aw = Op(a). 
Weyl was led to this choice because of the immediate ordering problems that 
occurred when one considered other observables than monomials a(x,p) = 
x k or a(x,p) = jr. For instance, using Schrodinger's rule what should the 
operator associated with a(x,p) = xp be? Weyl's rule immediately yields 
the symmetrized quantization rule 

a(X,P) = ±(XP + PX) 

and one finds that more generally (McCoy |27| . 1932) 

x m p n 1 (fj pn-k X mpk (4) 
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It turns out that the Weyl quantization rule for monomials is a 
particular case of the so-called "r-ordering" : for any real number r one 
defines 
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this rule reduces to Weyl's prescription when t = \. When r = 1 one 
gets the "normal ordering" X n P n familiar from the elementary theory of 
partial differential equations, and r = yields the "anti-normal ordering" 
p n X n sometimes used in physics. We now make the following fundamental 
observation: the Born-Jordan prescription ([I]) is obtained by averaging the 
r-ordering on the interval [0,1] (de Gosson [18], de Gosson and Luef [2T]). 
In fact 

k\(n-k)\ 
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One interesting feature of the quantization rules above is the following: 
suppose that the operators X and P are such that 

[x, P] = xp - Px = in. 

then [X m ,P n ] is independent of the choice of quantization; in fact (see 
Crehan [7] and the references therein): 
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[x m ,p n }= >^ (^) fe ( 7) [ n \p n - k x m - k . (7) 



In physics as well as in mathematics, the question of a "good" choice 
of quantization is more than just academic. For instance, different choices 
may lead to different spectral properties. The following example is due to 
Crehan [7]. Consider the Hamiltonian function 

H(z) = \{p 2 + x 2 ) + \{p 2 + x 2 f . 

The term that gives an ordering problem is evidently (p 2 +x 2 ) 3 ; Crehan then 
shows that the most general quantization invariant under the symplectic 
transformation (x,p) i — > (p, —x) is 

H = i(P 2 + X 2 ) + X(P 2 + X 2 ) 3 + X(3ah 2 - 4)(P 2 + X 2 ). 
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The eigenfunctions of H are those of the harmonic oscillator, and the cor- 
responding eigenvalues are the numbers 

E N = (N+ \)h + \h{2N + l) 3 + \h{2N + l)(3ah 2 - 4) 

(N = 0, 1, 2, ...) which clearly shows the dependence of the spectrum on the 
parameters a and A, and hence of the chosen quantization. 

3 Born— Jordan Quantization 

3.1 First definition 

Let Zq = (xq,Pq) and consider the "displacement" Hamiltonian function 
H zo = &(z,zo). The flow determined by the corresponding Hamilton equa- 
tions is given by f t (z) = z + tz ; for * G S'(IR 2n ) we define T(z )V(z) = 
(fi)*^>(z) = *&(z — zq). The r-quantization of H zo is the operator H ZQ = 
a(Z,zo), Z = (X,P); the solution of the corresponding Schrodinger equa- 
tion at time t = 1 with initial condition tp is given by the Heisenberg operator 
f(z ) = eTi^^o). its ac tion on ^ G S'(R n ) is explicitly given by 

f (z )i/}(x) = ei^-^^iPix - x ). (8) 

Let a G 5'(IR 2n ) be an observable (or "symbol"). By definition, the Weyl 
correspondence a Ay? is defined by 

%^ = (2^)" / a a (z)f(z^dz (9) 
J1 2 " 

where a a = F a a is the symplectic Fourier transform of a, that is 

for a G <S(M 2 ™); informally 

= (kr)" / e-i^a(z')dz'. (11) 

The symplectic Fourier transform T a is an involution (J- -2 = Jj); it is related 
to the usual Fourier transform T on M. 2n by the formula T a a(z) = Ta(Jz), 
J the standard symplectic matrix. The action of the operator on a 
function ip G <S'(IR n ) is given by 

A ^ = i^nT I a a (z )f(z )^dz . (12) 

JR 2n 
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How can we modify this formula to define Born- Jordan quantization? An 
apparently easy answer would be to first define r-quantization by replacing 
H ZQ by its r-quantized version H ZOjT , and then to average the associated 
operators T t (zq) thus obtained to get a "Tbj(zo) operator" which would 
allow to define j4bj- However, such a procedure trivially fails, because all 
r-quantizations of the displacement Hamiltonian H Zo coincide with H ZQ as 
can be verified using the polynomial rule ([5|). There is however a simple 
way out of this difficulty; it consists in replacing, as we did in [21J, T(zq) by 
Q(zo)T(zo) where 

= Bin ^°/f) (13) 
and we define the Born- Jordan operator Abj by 

= (2^)" / a a (z)e(z)T(z )i>dz. (14) 

jR 2n 

This formula will be justified below. 

3.2 Pseudo-differential formulation 

There is another way to describe Born-Jordan quantization. Writing for- 
mula (I12p in pseudo-differential form yields the usual formal expression 

A ^( x ) = / e^-^a(i(x + y),p)^y)dpdy (15) 

for the Weyl correspondence (we assume for simplicity that a € 5(M 2ra ) and 
ip G <S(R n )). We now define the r-dependent operator a la Shubin [30] : 



= {^rf I e^ x -^a(rx + (1 - r)y),p)^{y)dpdy- 

JR2n 



(16) 



the Born-Jordan operator A^j with symbol a is then defined by the average 

A B3 ^ = [ A T ipdT (17) 
J 

which we can write, interchanging the order of the integrations, 

Abj1>(x) = (2^)™ / eTM*-y) aB j(x,y,p)ij(y)dpdy (18) 



where 

aBj(x,y,p)= a(rx + (l-r)y,p)dT. (19) 
io 
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We have been a little bit sloppy in writing the (usually divergent) integrals 
above, but all three definitions become rigorous if we view the operators 
Aw, A T , and A B j as being defined by the distributional kernels 

K(x,y) = {^nf'i^a^ + ylP) (20) 

Krfrv) = (^l)" 72 PT^XC™ + 1 - r)y,p) (21) 

(J 7 ^ 1 is the inverse partial Fourier transform with respect to the second set 
of variables) and 

K B j(x,y) = [ K T (x,y)dT (22) 
Jo 

where J 7 ^ 1 is the inverse Fourier transform in the second set of variables. 
We will give below an alternative rigorous definition, but let us first check 
that definition (|18j ) (|22[) coincides with the one given in previous subsection. 
Define the modified Heisenberg-Weyl operators 

f T (z )4>(x) = e^ 2T - 1)poX °f(z )^(x) (23) 

that is 

T t (zq)iI>(x) = eK^-a-rJpoxo)^,^ _ XQ ). (24) 
These obey the same commutation rules 

T T (zo)T T (zi) = e^^f T ( Zl )f T (z ) (25) 
as the usual Heisenberg operators T(zq). 

Proposition 2 Let a G S'(M. 2n ), -0 £ <S(K n ). The Born-Jordan operator 
\17\) is given by formula {1$ , that is 

A B34> = (wiY / a a (z)f BJ (z)i;dz (26) 

with 

f w(2) _ eW f W , e W - ^. (27) 

In particular, A B j is the Weyl operator with symbol 

«BJ= (^l) n a*-^e (28) 
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Proof. (Cf. [21\ f!8j). One verifies by a straightforward computation that 
the Shubin formula (|16p can be rewritten as 

A T tp = / a a (z)f T (z)ipdz. (29) 

Let us now average in r over the interval [0,1]; interchanging the order of 
integrations and using the trivial identity 

/^-^r = ^sin^ 
Jo px 2h 

we get 

f BJ (z)= I T T (z)dr = @{z)f (z) 
Jo 

hence formula (126|) . To prove the last statement we note that formula (126ft 
can be rewritten 

= {^nT [ a (J (z)Q{z)f{z)i,dz 

= (sk)" / (aBj) CT (2)f(z)^dz 
in 2 " 

where (a-Qj)a = a©- Taking the inverse Fourier transform we get, noting 
that Ta{a *b) = (27rh) n J 7 rT aJ r r jb, 

hence ([28]) since = F a ® because the function G is even. ■ 

One easily verifies that the (formal) adjoint of A T = Op T (a) is given by 

Op T (a)* = Pl _ T (a) (30) 

and hence 

OpBj(a)* = Op BJ (a). (31) 

Born-Jordan operators thus share with Weyl operators the property of being 
(essentially) self-adjoint if and only if their symbol is real. This property 
makes Born-Jordan prescription a good candidate for physical quantization, 
while Shubin quantization should be rejected being unphysical for t^. 
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4 The Born— Jordan— Wigner Distribution 

4.1 The r- Wigner distribution 

In a recent series of papers Boggiatto and his collaborators [H EJ [3] have 
introduced a r-dependent Wigner distribution Wig r (/, g) which they aver- 
age over the values of r in the interval [0,1]. This procedure leads to an 
element of the Cohen class [22], i.e. to a transform of the type C(ip,(f)) = 
Wig T (ip, (f))*9 where 6 £ S'(R 2n ). From the point of view of time-frequency 
analysis this can be interpreted as the application of a filter to the Wigner 
transform. 

Let us define the r- Wigner cross-distribution Wig T (tp, (ft) of a pair (ip, 
of functions in iS(R n ): 

Wig T (^, </>)(z) = (^) " [ e-TiPV^x + Ty)4>(x - (1 - r)y)dy. (32) 

Choosing t = \ one recovers the usual cross- Wigner transform 

Wig(V, 4>){z) = (^i)" / e-^Mx + \y${x - \y)dy (33) 
and when r = we get the Rihaczek-Kirkwood distribution 

fl(V,0(z) = {^n) n ' 2 e>^(p) (34) 

well-known from time- frequency analysis [22] . 

The mapping Wig r is a bilinear and continuous mapping 5(IR n ) x5(M n ) — 
<S(M 2n ). When ip = <fi one writes Wig^^VO = Wig T ^; it is the r- Wigner 
distribution considered by Boggiatto et al. [H El [3]). It follows from the 
definition of Wig T that we have 

Wig r (<^) =Wi gl _ T (^); (35) 

in particular 

Wig7? = Wig 1 _ r V (36) 
hence Wig r ip is not a real function in general if r ^ ^ . 

Proposition 3 Assume that ip,(/> E L 1 (R n ) n L 2 (R n ) . Then 

( Wg T (if> ) </>)(z)dp = if>(x)^(x) (37) 
it" 

and 

[ Wig T ^,<j>)(z)dx = T^{p)F4>{p). (38) 
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Proof. Formula (j37[) is straightforward. On the other hand 

Wig T (^, 4>){z)dx = (2^)" / e~^ py i(;(x + ry)4>(x - (1 - r)y)dxdy 
and setting x' = x + Ty, x" = x — (1 — r)y we have dx'dx" = dxdy so that 
Wig T (V, 0)(z)cfc = (^) n / e'y x '^{x')e^{x")dxdy 



hence formula (|38|) . Notice that the right-hand sides of (|37j) and (|38|) are 
independent of the parameter r. ■ 

In particular [T], the r-Wigner distribution Wig r ip = Wig r (V>,V0 satis- 
fies the usual marginal properties: 

Wig T ^{z)dp = \^(x)\ 2 , [ Wig T ^{z)dx = \Fi){p)\ 2 . (39) 

There is a fundamental relation between Weyl pseudo-differential oper- 
ators and the cross- Wigner transform, that relation is often used to define 
the Weyl operator Ay] = Op w (a): 

<AvV#) = (a,Wig(iM)) (40) 

for if),(f) € 5(M ra ). Not very surprisingly this formula extends to the case of 
r-operators: 

Proposition 4 Let ip,<ft € 5(M n ), a G 5(]R 2n ), and r a real number. We 
have 

(A r ^) = (o,Wig T (^,0)) (41) 
where (•, •) is i/ie distributional bracket on M 2n and A T = Op T (a). 

Proof. By definition of Wig r we have 
{<i,Wig r = 

Defining new variables x' = x — (1 — r)y and y' = x + ry we have y = y' — x' , 
dydx = dy'dx' and hence 

<a,Wig T (V^)) = 

(lil)" /. e-^'^')a(rx'+(l-r)n',p)^(n')0(x / )<iy / a'p'a'x / ; 
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the equality f)41[) follows in view of definition (|16p of A T . ■ 

Formula (|4ip allows us to define = Op T (a)ip for arbitrary symbols 
a € 5'(lR 2n ) and ^ € 5(lR n ) in the same way as is done for Weyl pseudo- 
differential operators: choose <fi G <S(M n ); then Wig T (ip, </>) £ <S(M 2n ) and the 
distributional bracket (a, Wig T (^>, </>)) is thus well-defined. This defines A T 
as a continuous operator tS(M n ) — > S'iW 1 ). 

4.2 Averaging over r 

We define the (cross) Born-Jordan-Wigner (BJW) distribution of ip, (j) € 
5(M n ) by the formula 

Wig B j(^0)(z)= / Wig T (V,^r. (42) 
J o 

We set Wig BJ ^ = Wig BJ (V>, V0- The properties of the BJW distribution are 
readily deduced from those of the r-Wigner distribution studied above. In 
particular, the marginal properties (137p and (138p are obviously preserved: 

/ Wig BJ (^0)(*)dp = ^(a;)0(aO; (43) 

and 

/ Wig BJ (^)(*)cfe = Tif>(pjjy@j. (44) 

JR™ 

An object closely related to the (cross-)Wigner distribution is the (cross- 
ambiguity function of a pair of functions ^,<j> £ <S(M n ): 

Aty, = (^-)" / e-^'^x' + - \x)dx'. 

JR« 

It turns out that A(tp, </>) and Wig T (^, 4>) are obtained from each other by a 
symplectic Fourier transform: 

Amb(V>, 4>) = T a Wig(V>, 4>) , Wig(V>, 4>) = T a Amb(V>, </>) 

{J- a is an involution), thus justifying the following definition: 

Amb BJ (V>, <j)) = T a Amb(V>, </>)• 

An important property is that the BJW distribution of a function ip is 
real, as is the usual Wigner distribution. In fact using the conjugacy formula 
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([36]) we have 



Wig BJ V = / Wig r ipdr = / Wig x _ T ipdr 
Jo Jo 
1 

Wig T ipdr = Wig BJ ip. 

Formula (|4ip relating Shubin's r-operators to the r- (cross) Wigner dis- 
tribution carries over to the Born- Jordan case: 

Corollary 5 Let a € S'(M. 2n ), A BJ = Op BJ (a). We have 

( J 4 BJ ^|0) = (a,Wig BJ (V,0)} (45) 

for all ip,(f>e S(R n ). 

Proof. It suffices to integrate the equality (A T ip\(j)) = (a, Wig T (?/>, </>)} with 
respect to r € [0, 1] and to use definitions (fT7|) and ([4*2"|) . ■ 

One defines Abj = Op BJ (a) for arbitrary a G 5'(IR 2n ) by the same 
procedure as for Weyl operators and noting that Wig BJ (?/>, 0) € 5(M 2n ) if 

V>,</> e <S(R n ). 

The following consequence of Proposition [4] will be essential in our study 
of the uncertainty principle: 

Proposition 6 Let be defined by (2l\ ). We have 

Wig BJ ty, </>) = Wigty, <j>) * TQ (46) 

and 

Amb BJ (^, <p) = {2Tih) n Amb(^, <j>)Q (47) 

Proof. In view of formula (|28p in Proposition [2] and formula (|42p above we 
have 

(o,WigBj(^,0)> = (6,Wig(^^)) 
where b = (27r/i) n a * J- a Q, hence 

(a, Wig B j(^, (/))) = (27r^f(a*7- (7 e,Wig^,0)} 

= (2nH) n (a, [T a Q o (-/ d )] * Wig(V>, </>)) 

(the last equality because is an even function). Since T a Q = TQ because 
G is even and invariant under permutation of the x and p variables, this 
proves formula ([46 p . and formula (|47p follows, taking the symplectic Fourier 
transform of both sides. ■ 
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4.3 Symplectic (non-)covariance 

he symplectic group Sp(2n, M.) is by definition the group of all linear auto- 
morphisms s of R 2n which preserve the symplectic form a(z, z') = Jz ■ z'\ 
equivalently s T Js = J. The group Sp(2n,R) is a connected Lie group, and 
its double covering Sp 2 (2n,R) has a faithful (but reducible) representation 
by a group of unitary operator, the metaplectic group Mp(2n,R) (see Fol- 
land [10], de Gosson [El [16]). That group is generated by the operators J, 
M L , m , and V- P defined by J = e ~ in ^^T and 

M L , m ^{x) = i m y/\detL\i/;(Lx) , % = e~m Px - x ip( x ) 

where L £ G£(n,M) and P E Sym(n,M); the integer m corresponds to the 
choice of an argument of det L. Denoting by 7r Mp the covering projection 
Mp(2n,M) — > Sp(2n,M) we have tt m p(J) = J and 

vr Mp (M,, m ) = (V L ° T ) , ^(Vp) =(y J) . 

Let now Aw = Op w (a) be an arbitrary Weyl operator, and s £ Sp(2n, R). 
We have 

SAwS' 1 = Op w (a o s- 1 ) (48) 

where S € Mp(2n,M) is anyone of the two metaplectic operators such that 
7r Mp (S) = s. This property is really characteristic of the Weyl correspon- 
dence; it is proven [TU [16] using the identity 

ST(z) = f(sz)S (49) 

where T[z) is the Heisenberg operator. One can show that if a pseudo- 
differential correspondence a < — > Op(a) (admissible in the sense above, 
or not) is such that Op(a o s _1 ) = S Op(a)S^ 1 then it must be the Weyl 
correspondence. For Born- Jordan correspondence we do still have a residual 
symplectic covariance, namely: 

Proposition 7 Let Abj = Op BJ (a) with a € S'(R 2n ). We have 

50p BJ (a) < S- 1 = Op BJ (ao S - 1 ) (50) 
for every S in the subgroup o/Mp(2n,R) generated by the operators J and 

M Lt m (With 7T M P(S) =S). 
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Proof. (Cf. de Gosson [IB]). It suffices to prove formula ([50]) for S = J 

and S = Mi^ m . Let first S be an arbitrary element of Mp(2n,M); we have 



SOp BJ (a) = (^) n J a a {z)Q(z)ST(z)dz 

= (2^)" (/ a a (z)@(z)T(sz)dz} S 

where the second equality follows from the symplectic covariance property 
(|49p of the Heisenberg operators. Making the change of variables z' = sz in 
the integral we get, since det s = 1, 



a a {z)Q{z)T{sz)dz = J a^s" 1 z)6(s" 1 z)T(z)dz. 
Now, by definition of the symplectic Fourier transform we have 

^{s~ l z) = (^) n J e-rA s ' lz ^a(z')dz' = (a o 8 -\(z). 

Let now S = Mr, m ; we have 

similarly G(J _1 z) = Q(z), hence in both cases 

SO PBJ (a) = {^- h ) n (| (a o s -%e(z)T(z)dz\ S 
= OpBj(aos _1 )5 
whence formula (1501). ■ 



5 The Uncertainty Principle 

We begin by reviewing the notion of density matrix (or operator) familiar 
from statistical quantum mechanics. The notion goes back to John von 
Neumann [32] in 1927, and is intimately related to the notion of mixed state 
(whose study mathematically belongs to the theory of C*-algebras via the 
GNS construction). This will provide us with all the necessary tools for 
comparing the uncertainty relations in the Weyl and Born-Jordan case. 
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5.1 Density matrices 

A density matrix on a Hilbert space H is a self-adjoint positive operator on 
T~L with trace one. In particular, it is a compact operator. Physically density 
matrices represent statistical mixtures of pure states, as explicitly detailed 
below. 

We will need the two following results: 

Lemma 8 A self- adjoint trace class operator p on a Hilbert space H is 
a density matrix if and only if there exists a sequence (ay)j>i of positive 
numbers and a sequence of pairwise orthogonal finite- dimensional subspaces 
{ / Hj)j>i of H such that 

where Hj is the orthogonal projection of % on 7~Lj and rrij = dim Hj 

It is a consequence of the spectral decomposition theorem for compact 
operators (for a detailed proof see e.g. [TE], §13.1). 

Lemma 9 Letip G L 2 (R n ), ip / 0. The projection operator P^ : L 2 (IR n ) — ► 
: A € C} has Weyl and Born-Jordan symbols given by, respectively 

p w = (27r^)"Wig w V (51) 

and 

p B j = (2^rWig BJ V>. (52) 
In particular and pbj o,re real functions. 

Proof. We have P^cj) = (<?!>|V')L 2 ' i / ; hence the kernel of Pw, is = ip <8> ip. 
Using a Fourier transform formula (|2ip implies that the r-symbol p T of P^ 
is given by 



Pr(z) = I e hPVK^{x + ry,x - (1 - r)y)dy 
e~h pv ip(x + ry)ip(x — (1 — r)y)dy 

= (2^/i) n Wig r V^)- 

Setting r = \ we get formula ()5ip . Formula (|52p is obtained by integrating 
p T (z) with respect to r € [0, 1]. That pw and pbj are real follows from the 
fact that both the Wigner and the WBJ distribution are real. ■ 

The following result describes density matrices in both the Weyl and 
Born-Jordan case in terms of the Wigner formalism: 
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Proposition 10 Let p be a density matrix on L 2 (M n ). There exists an or- 
thonormal system (ipj)j>i of L 2 (M. n ) and a sequence of non-negative numbers 
(Aj)j>i such that Y2j>i -\j = 1 an d 

P = °Pbj Pbj = Op w pyj (53) 

the symbols p\y and pbj being given by 

p w = (27rR) n ^A i WigVi (54) 
i>i 

and 

/ 9Bj = (2vr^r^A J Wig B j^- (55) 

Proof. Taking % = L 2 (IR n ) in Lemma [5] we can write p = ^ - a^IL; 
where each ILy is the projection operator on a finite dimensional space 
Hj C L 2 (IR n ), and two spaces Uj and Hi are orthonormal if j ^ £. For 
each index j let us choose an orthonormal basis Bj = (ipj+i, ipj+m } ) of 
"Hj ; the union B = DjBj is then an orthonormal basis of (BjHj, and we 
have, using Lemma [8] 

j>l \j+l<k<j+mj j 

where EL^ is the orthogonal projection on the ray {Xijjj : A € C}. Since 
each index ay is repeated mj times due to the expression between brackets, 
this can be rewritten 

with the Xj = mjCtj summing up to one. In view of Lemma[9]the Weyl (resp. 
Born-Jordan) symbol of IL^ is (2-Kh) n Wig w tp (resp. (2irh) n Wig BJ i/j) 
hence the result. ■ 

Notice that the orthonormal bases Bj in the proof can be chosen arbitrar- 
ily; the decompositions (|54p and ([55]) are therefore not unique. (In Physics, 
one would say that a mixed quantum state can be written in infinitely many 
way as a superposition of pure states, a pure state being a density operator 
with symbol a Wigner function). 
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5.2 A general uncertainty principle 

In what follows the notation 

a < — ► A = Op(a) 

is indifferently the Weyl or the Born- Jordan correspondence. Both have the 
property: 

If the symbol a is real, then the operator A is essentially self- 
adjoint (in which case we call it an observable). 

Notice that the Shubin correspondence does not have this property for 
t^j since A* = A\_ T . 

Let be a density matrix on L 2 (R n ). We assume that p = (2irh) n Op(p). 
The function p is a real function on phase space M? n and we have 

Trp = / p(z)dz = 1. (56) 

Observe that we do not in general have p > 0. 

Let A be an observable. Its expectation value with respect to p is by 
definition the real number 

(A) = [ a(z)p(z)dz (57) 

where it is assumed that the integral on the right side is absolutely conver- 
gent. We will write, with some abuse of notation, 

(A) = Tv(pA) (58) 

(see the discussion in de Gosson |16j . §12.3, of the validity of various "trace 
formulas" ) . In view of formulas (|40p and (|45p we have either 

(l) = (2 7 r^) n ^A i (a,Wig^) 

(when a < — > A is the Weyl correspondence) or 

(l) = (2^r^A J (a,Wig BJ ^> 
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(when a < — > A is the Born-Jordan correspondence) and hence 

(A) = (2^)" £ XiiMMAv = *M)j (59) 

where (A)j = (Aipj\ipj) L 2 (recall that {tpj)j>i is an orthonormal system). 
If A 2 also is an observable and if (A 2 ) = Tv(pA 2 ) exists, then the number 

(V a r p A) 2 = (A 2 ) p -(A) 2 p (60) 

that is 

(Var p A) 2 = Ti(pA 2 ) - Tv(pA) 2 (61) 

is the variance of A; its positive square root Var p A is called "standard devi- 
ation". More generally consider a second observable B; then the covariance 
of the pair (A, B) with respect to p is defined by 

Cov p (A, B) = Tr(pAB) - Tr(pA) Ti(pB). (62) 

It is in general a complex number, and we have 

Cav p {A, B) = Cov p (B, A). (63) 

The covariance has the properties of a complex scalar product; it therefore 
satisfies the Cauchy-Schwarz inequality 

Ti(pA 2 ) Tr(pB 2 ) > | Cov p (A, B)\ 2 . (64) 

The following lemma will be useful in the proof of the uncertainty in- 
equalities below: 

Lemma 11 If the covariance of two observables A and B exist we have 

Re Cov p (A, B) = \ Tr(p{A, B}) - Tr(pl) Tr(^B) (65) 
1mCav p (A,B) = j i Tr(p[A,B]) (66) 

where {A, B} = AB + BA and [A, B] = AB — BA are, respectively, the 
anticommutator and the commutator of A and B. 
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Proof. We have, in view of (|63|). 

2 Re Cav p (A, B) = Cov p (A, B) + Cov p (£, A) 

= Tr(p{A, B}) - 2 Tr(pl) Tr(pS); 

for the second equality 

2z Im Cov p (l, B) = Cov p (l, 5) - Cov p (i?, 1) 
= Tr(pAB) - Tr (pB A) 
= Tr(p[A,B}). 

u 

Observe that the anticommutator and commutator obey the relations 

{A,B}* = {A,B} , [A, B]* = —[A, B] 

and hence ([A, B]) is a pure imaginary number or zero; in particular \([A, B])\ 2 < 
0. 

Proposition 12 // the variances and covariance of two observables A and 
B exist then: 

(Var p A) 2 (Var p 5) 2 > CovJ m (A,B) 2 - \{[A,B]) 2 p (67) 
where ([A,B]) 2 < and 

CovJ^iA, B) = ±(Cov„(A, B) + Cov p (£, A)) (68) 
is a real number. In particular the Heisenberg inequality 

(Var p A) 2 (Var p 5) 2 >4([l,^ (69) 

holds. 

Proof. Replacing A and B with A — {A) and B — (B) it is sufficient to 
prove ([6?]) when (A) = (B) = 0. We thus have to prove the inequality 

TripA 2 ) Tr(pB 2 ) > Cov p ym (l, B) 2 - \ Tr(pL4, B}) 2 . (70) 

Noting that definition (|68p can be rewritten 

Cov^ m (A, B) = \ Tr(p{l, B}) - Tv(pA) Tv(pB). 
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we remark that in view of formulas (|65p and (|66p in the lemma above we 
have 

| Cor p (A, B)\ 2 = \ Tr(p{A, B}) 2 - \ Tr(pL4, B]f (71) 
= Cov7 m (l, Bf - \ Dr(p[A, B]) 2 (72) 

hence the proof of (|7U|) is reduced to the proof of the inequality 

TrG?! 2 ) Tr(pB 2 ) > \Cov p {A,B)\ 2 (73) 

which is just the Cauchy-Schwarz inequality (f64"P for covariances. ■ 

5.3 Weyl vs Born— Jordan 

Let us discuss the similarities and differences between the uncertainty prin- 
ciples associated with the Weyl and Born-Jordan correspondences. First, 
as already observed in the Introduction, the Weyl and Born-Jordan quan- 
tizations of monomials x^p" (and hence of their linear combinations) are 
identical when m + n < 2. This implies, in particular, that if the symbols 
a and b are, respectively, multiplication by the coordinates Xj and pj then 
the corresponding operators A and B are, in both cases given by Xj = Xj 
and Pj = —ihd Xj . It follows that V&r p Xj and Var p Pj satisfy the usual 
Robertson-Schrodinger inequalities 

Var p Xj Var p Pj > CovJ m (Xj, Pj) 2 + \h 2 . (74) 

We mention that the inequalities (|74p can be rewritten in compact form 

as 

S + \ihJ > (75) 

where > means "semi-definite positive", J is the standard symplectic 
matrix, and 

(Cov p (X,X) Cov p (X,P)\ 

\Cov p (P,X) Cov p (P,P)J { 1 

is the statistical covariance matrix. The formulation (|75p of the Robertson- 
Schrodinger inequalities clearly shows one of the main features, namely the 
symplectic covariance of these inequalities, which we have used in previous 
work [171 [20] to express the uncertainty principle in terms of the notion 
of symplectic capacity, which is closely related to Gromov's non-squeezing 
theorem from symplectic topology. This has also given us the opportunity 
to discuss the relations between classical and quantum mechanics in [19J. 
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Let S € Mp(2n,R) i S = tt m p(S) and set A' = SAS' 1 , B' = SBS' 1 ; we 
are assuming that A,B correspond, as in the proof of Proposition [T2"l to an 
arbitrary quantization scheme a i — > A. We have quite generally, using the 
cyclicity of the trace, 

Cov^ m (A, B) = \ Tr(p{A, B}) - Tr(pA) Ti(pB) 

= \ TiipS^iA', B'}S) - TrCpS- 1 A' S) TrtfS^B'S) 
= \ TriS^pSiA', B'}) - TviS^pSA') TriSpS^B') 
= Cov|- M',B>); 

similarly Var p A = Var^_ x A and {[A,B]) 2 p = ([A', B'})~^~ _ 1 and hence 

(Va V A') 2 (Va V B') 2 > CovJT^, B'f - ±<[A', B%. 

with p' = SpS -1 . Suppose now that the operator correspondence a < — > A 
is the Weyl correspondence; then, by Proposition [7] we have SpS^ 1 = Op(po 
s _1 ) and the inequalities ([67]) become 

(Var^-x A^^Var^-! £?W) 2 > Cov^AW,^) 2 " K^w.^wD^-i- 

(77) 

Again, in view of Proposition [7J in the Born-Jordan case we have inequality 
(Var^A'ej^Var^ 

~ (78) 
only for those S € Mp(2n,R) which are products of metaplectic operators 
of the type J and Ml jTu . 

6 Discussion 

There is an old ongoing debate in quantum mechanics on which quantization 
scheme is the most adequate for physical applications; an interesting recent 
contribution is that of Kauffmann |25| , who seems to favor the Born- Jordan 
correspondence. The introduction of the r-Wigner and Born-Jordan distri- 
butions has been motivated in time-frequency analysis by the fact that the 
usual cross- Wigner distribution gives raise to disturbing ghost frequencies; 
it was discovered by Boggiatto and his collaborators [DEI [3] that these ghost 
frequencies were attenuated by averaging over r. 

The study of uncertainties for non-standard situations has been tackled 
(from a very different point of view) by Korn [26]; also see the review pa- 
per [TT] by Folland and Sitaram, which however unfortunately deliberately 
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ignores the fundamental issue of covariance. Gibilisco and his collaborators 
|12l [T3l Q3] give highly nontrivial refinements of uncertainty relations us- 
ing convexity properties, and studied the notion of statistical covariance in 
depth. 

We mention that in a very well written thesis, published as a book, 
Steiger [31] has given an interesting historical review and analysis of the 
evolution of the uncertainty principle; in addition he compares the interest 
of several different formulations, and gives a clever elementary derivation of 
the Robertson-Schrodinger inequalities for operators. The work also con- 
tains Matematica codes for the computation of (co-)variances. 

Acknowledgement 13 This work has been supported by a research grant 
from the Austrian Research Agency FWF (Projektnummer P23902-N13). 
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